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1.  INTRODUCTION 

Under  the  assumption  of  high  frequency,  wave  propagation  can  be  described  by 
the  ray  method.  Many  quantities  in  wave  phenomena  are  characterized  by  some 
ray  parameters  and  their  derivatives.  For  example,  the  slope  of  arrival  time-offset 
curve  is  characterized  by  the  emergence  angle  of  a  ray;  the  curvature  is  characterized 
by  derivative  of  that  angle,  and  the  amplitude  of  the  signal  is  characterized  by  the 
derivative  of  the  taking-off  angle.  Moreover,  when  we  solve  an  equation  containing 
ray  parameter,  we  also  require  analytic  expressions  for  the  derivatives.  Generally, 
it  is  difficult  to  obtain  analytic  expressions  for  these  derivatives.  Wc  can  use  finite 
difference  approximation  for  these  derivatives.  This  type  of  approximation  fails  at 
caustics.  And  also  the  accuracy  and  stability  of  the  calculation  are  influenced  by  the 
smoothness  of  interfaces  and  the  spacing  between  adjacent  rays. 

In  this  paper,  we  give  a  recursive  algorithm  for  ray-parameter  derivatives,  when 
the  medium  is  made  up  of  constant-velocity  layers  separated  by  arbitrary,  smooth 
interfaces.  The  algorithm  uses  derivatives  of  interfaces,  intersection  points  of  rays 
with  the  interfaces,  and  angles  of  incidence  and  transmission,  all  of  which  are  easily 
computed.  A  conclusion  of  these  formulas  is  that  no  ray  caustic  exists  for  transmission 
or  for  the  reflection  from  any  one  of  the  reflectors  when  the  interfaces  are  flat. 

2.  RECURSIVE  ALGORITHM 

Suppose  that  an  acoustic  earth  model  is  made  up  of  constant-velocity  layers  sep¬ 
arated  by  arbitrary,  smooth  interfaces.  We  are  going  to  derive  a  recursive  formula  for 
computing  ray-parameter  derivatives  when  a  ray  propagates  through  the  interfaces. 

We  consider  two  interfaces  L\  and  £2,  defined  respectively  by  r  =  /j(.r)  and 
z  =  /2(i).  A  ray  from  L\  arrives  at  L2,  with  intersection  points  (xi,ri)  and  ro), 
The  incident  angle  from  the  vertical  at  L\  is  cq  and  the  outgoing  angle  is  o_>.  The 
angle  of  the  normal  at  L\  is  6.  When  the  ray  arrives  to  L 1,  there  two  possibilities: 
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reflection  and  transmission.  Suppose  that  for  this  ray  the  initial  angle  from  the 
vertical  is  P  and  the  initial  depth  is  r0. 

Starting  from  the  tangent  equation 

~lx  lr=I1=  tan6>’ 


the  Snell'  law 


sin(ai  —  9)  sin(a2  -  9) 


for  transmission, 


Ci  c2 

a2  -  9  =  9  -  for  reflection, 

and  the  stright  line  equation 

x2-  xi  =  (z2  -  z\ )  tan  a2, 

and  differentiating  them  with  respect  to  /?,  we  find  a  recursive  algorithm  as  follows: 


or 


d9  cPfi  dxi2 

M  =  -l?l“’'70cos  6 ■ 

da2  d9  c2cos(oi  —  9)  ,dot\  d9 

——  =  —  + - t - —  (— -r  — — )  for  transmission, 

dp  dp  ci  cos(a2  -  9y  dp  dp 

da2  nd9  doi\  , 

— —  =  2  — - —  for  reflection, 

dp  dp  dp 


dfr 

dx 


dx2  z2  —  2i  da2 


dfi 


dx  i 


(1  _  tan ot2—  \X=I2)—  =  77^“—  +  U  ~  tana2— 

'■‘t* 


dP  cos2  a2  dP 

The  initial  values  of  the  relevant  parameters  are 

da  i  _ 
dp  ~  ' 

dx i  Z\  —  zq  f ,  df i  *  . 

_  =  — - — /(l  -  tan  Qi  —  I=XI). 
dp  cos2  ai  dx 

3.  APPLICATION 


(1) 

(2) 

(3) 

(4) 


Geometry  spreading  factor  OP  /  d£ 

From  Docherty  (1987),  we  require  the  calculation  of  dp  /  d£  in  order  to  obtain 
the  amplitude  of  the  signal.  Here, 

dp_  _  .dx 2  , 
di  ~  • ' dp  113=0 
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where  3  is  the  taking-off  angle  of  the  ray  and  <f  is  the  horizontal  coordinate  of  the 
end  of  the  ray.  The  factor  dx^/dp  can  be  calculated  by  using  equation  (4). 


Curvature  of  arrival  time-offset  curve 


Starting  from  the  equation 

dr  sin  a 
dx  v 

where  a  is  the  emergence  angle  of  the  ray  and  v  is  the  velocity,  it  follows  that 

d2r  cos  a  da 
Ox2  v  dx 


Since  the  upper  surface  is  horizontal,  a  =  c*2  at  z  =  0.  From  equations  (2)  to  (4),  we 
can  compute  da^ldP  and  dxo/dP  from  which  we  find 


da 

dx 


da?  .  .  dx-2  , 

dp  |q2=q  '  dp  Ui=r 


Caustic  for  flat  interfaces 

An  interesting  conclusion  of  equations  (1)  to  (4)  is  that  when  interfaces  arc  Hat. 
no  caustic  exists,  i.e.,  A P  ^  0  implies  Ax  ^  0.  In  fact,  for  flat  interfaces,  dQ/dp  =  0 
from  equation  (1).  For  the  transmission,  all  coefficients  are  positive  in  equations  (2) 
and  (4),  so  dx-ijdP  is  positive.  For  the  reflection,  dp 2  and  z2  -  zq  change  the  sign  at 
same  time,  so  cfx2/d/?  still  keeps  the  positive  sign.  This  implies  that  rays  from  two 
different  taking-off  angles  must  not  intersect  each  other  for  transmission  or  leflection 
from  one  reflector. 


4.  CONCLUSION 

The  recursive  algorithm  has  higher  accuracy  and  more  stability  than  the  finite 
difference  approximation.  Compared  the  ray  tracing,  only  a  little  additional  compu¬ 
tation  is  required  for  the  recursive  algorithm.  Furthermore,  because  of  its  analytic 
expression,  the  recursive  algorithm  provides  a  useful  tool  for  the  theoretical  analysis. 
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Transmission 


Fig.  1.  Paths  of  arriving  ray  for  transmission  and  reflection. 
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